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Abstract 

Cosmic ray nuclei, cosmic ray electrons with energy above a few GeV, and the 
diffuse gamma-ray background radiation (GBR) above a few MeV, presumed to be 
extragalactic, could all have their origin or residence in our galaxy and its halo. The 
mechanism accelerating hadrons and electrons is the same, the electron spectrum 
is modulated by inverse Compton scattering on starlight and on the microwave 
background radiation; the 7-rays are the resulting recoiling photons. The spectral 
indices of the cosmic-ray electrons and of the GBR, calculated on this simple basis, 
agree with observations. The angular dependence and the approximate magnitude 
of the GBR are also explained. 



PACS numbers: 98.70.Sa, 98.70.Rz, 98.70.Vc. 



The origin and spectrum of non-solar cosmic ray hadrons, electrons and photons 
has been debated for almost a century. The cosmic ray (CR) nuclei have a power-law 
spectral flux dF/dE oc E~ ft with a series of break-point energies (/?i,2,3 ~ 2.7, 3.0, 
and 2.5 in the intervals 10 10 eV < E < E knco ~ 3 x 10 15 eV, E knee < E < E ankle 
~ 3 x 10 18 eV, and E ankle < E < E GZK ~ 6 x 10 19 eV). Below E knee , protons constitute 
~ 96% of the nucleon flux at fixed energy per nucleon. Their flux is [|TJ: 

-2.70±0.05 

cm" 2 s" 1 sr" 1 GeV _1 . (1) 

A presumed extragalactic origin |2j of the highest-energy CRs is at variance with 
observations |§ of CRs above the 'GZK cutoff': 7-rays and neutrinos are not good 
candidate primaries J| and, for E > Egzk, the flux of protons and nuclei originating 
more than ~ 20 Mpc away || should be quenched || by the cosmic background 
radiation (CBR). In this letter we argue that, above typical nuclear energies of a few 
MeV, non-solar CR electrons and the 'diffuse' 7-ray background may be intimately 
related to hadronic CRs, whose observed spectral shape we use as input. 

By assuming a common accelerating mechanism for all charged particles we de- 
duce a power law that fits, above E e = 5 GeV, the CR electron flux |7J || shown in 
Fig. 1, from the hadron spectral index j3\. Unlike protons, electrons are significantly 
affected by Compton scattering on starlight and on the CBR. The up-scattered pho- 
tons have a spectral index that agrees with that of the observed [|[ diffuse 7-ray 
background above a few MeV, shown in Fig. 2. To obtain the observed 7-ray inten- 
sity, we need a scale-height of the CR electrons much bigger than the accepted ~ 1 
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kpc minimum |10| . There is evidence for CR electrons well above galactic disks |l I 
from observations of synchrotron emission by galaxies seen edge-on. We envisage 
the possibility that the 'box' in which charged CRs reside be much larger than the 
galactic disk, and we call this box 'the halo' [ |12|j . 

A classic argument determines the volume of a region to which CRs are confined. 
An abundance ratio of secondary to primary rays (e.g., Li/C) is used to deduce X: 
the mean column density they traverse. Ratios of unstable to stable isotopes (e.g. 
10 Be/ 9 Be or 26 A1/ 27 A1) are used to deduce the mean length, L, of their voyage. The 
inferred mean density [|13| is p = X/L ~ p(disk)/5. The usual conclusion -that the 
confining volume is ~ 5 V(disk)- relies on tacit assumptions that need not be right. 
If the halo is magnetized, as the disk is, we have a "leaky box" within a (much less 
dense) leaky box. If the lifetimes of 10 Be and 26 Al are not longer than the mean 
residence time in the disk, the observed ratios determine L(disk), but not L(halo). 
The stable isotope ratios do determine the total column density X(tot), but this 
may include various prolonged journeys in the halo. Without knowing L(halo) the 
density p = X(tot)/L(tot) of the overall confining domain cannot be estimated. The 
conventional tacit assumption is that the stable CRs that reach us cannot have spent 
long periods in a baryonically thin halo, before or between disk-crossings. 

To relate the spectra of CR electrons and protons we shall need an estimate 
of the protons' spectrum at their source. At E < Ek nee a source spectrum dF s /dE 
with index j3 s ~ 2.2 is obtained from collisionless-shock simulations |]T3 or analytical 



estimates of acceleration by relativistic jets [HJ. The CR spectrum is modulated by 
their residence time in the galaxy, r gal (E). For a steady source of CRs the energy 
dependence of the observed flux is roughly that of r ga idF s /dE. Observations of 
astrophysical and solar plasmas and of nuclear abundances as functions of energy 
indicate flT§ that r gal (E) oc E -0 - 5 * 0,1 , explaining ft ~ 2.7, as in Eq. (Q). 

The CR electron flux of Fig. 1 is well fit, from ~ 10 GeV to ~ 2 TeV, by: 



,p r tti 1 -3.2±0.10 

e ~ (2.5 ±0.5) x 10 5 
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s" 1 sr" 1 MeV" 1 . (2) 



The e + admixture is ~ 7% from 3 to 50 GeV 17]. To a good approximation, as we 
assume at all energies, CR electrons are not pair produced, nor are they secondary 
products of interactions, nor the result of relic-particle annihilation or decay. 

Practically all CR acceleration mechanisms invoke a region of space that is swept 
by a moving magnetic field, such as would be carried by the rarefied plasma in a 



supernova shell @ or by a 'plasmoid' of jetted ejecta [12|. The magnetic field 
acts as a moving 'mirror' that imparts the same distribution in velocity, or Lorentz 
factor 7 = E/m c 2 , to all charged particles. To the extent that particle-specific losses 
(such as synchrotron radiation) can be neglected at the acceleration stage, all source 
fluxes have the same energy-dependence. Confinement effects preserve this equality 
for ultrarelativistic electrons and protons: their behaviour in a magnetic maze is 
the same. But, unlike for hadrons, the 'cooling' time of electrons in their 'inverse' 
Compton scattering (ICS) interactions with starlight and the CBR is shorter than 



their galactic confinement time, r ga i(E), above a relatively low energy. 

Consider electron interactions with the CBR. Let To = 2.728 K, no ~ 411 cm -3 , 
and €q w 2.7 kTo ~ 6.36 x 10~ 10 MeV be the current CBR temperature, number 
density and mean energy fI8| . The cross section for ICS is a T « 0.65 x 10~ 24 cm 2 
(for the relevant E e range the Thompson limit is accurate, even for ICS on starlight). 
The mean energy E 7 of the upscattered photons, -or AE e , the mean energy loss per 
collision- is E 7 m AE C « (4/3)eo(E e /m e c 2 ) 2 . The single electron interaction rate is 
cr T no c. Let R (also an inverse time) be the production rate of CR electrons, as- 
sumed to be constant JE], and let dn^/dEg be their source number- density spectrum. 
The actual density dn e /dE e in an interval dE e about E e is continuously replenished 
and depleted by electrons whose energy is being degraded by ICS. This leads to a 
steady-state situation in which production and losses are in balance: 

d / dn c \ ^ dn^ . . 

For a relatively uniform galactic CR electron density, Eq. (|3|) also applies to the 
local electron flux dF e ~ (c/47r)dn e . Substitute the spectrum dn|/dE ~ E"' 38 into 
the flux version of Eq. fl3|) to obtain: 

dF P 3m 2 c 4 R dF? 



e v " e 



cx E e "^ +1 ) . (4) 



dE c 4 (J3 B - 1) c <tt n e E e dE e 
Compton scattering on starlight and synchrotron radiation on magnetic fields can be 
included by replacing n e by the total energy density + B 2 / (Stt), whose local 
variations do not affect the spectral shape. For electrons with E e < (m e /m p ) Ek ne e~ 



1.6 TeV we deduced that (3 S ~ 2.2 . Thus, (3 S + 1 = 3.2, in agreement with the 
observed spectral index, Eq. (Q) and Fig. 1. Above the 'electron's knee' at E c ~ 1.6 



TeV the spectrum should steepen up by A/3 ~ 0.25, like that of CR hadrons |L5 
The available spectral measurements extend only to E c < 1.5 TeV. 



The cooling time of electrons in the CBR bath is given by 



fE x E c _ 3m e 2 c 4 

Ticsl e) ~ n a T cAE c " 4E e n eo<x 1 
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GeV 
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(5) 



At our location in the Galaxy, starlight and magnetic fields have energy densities 
similar to that of the CBR, and r ICS is shorter by a factor ~ 3. The galactic escape 
time of GeV electrons should be similar to that of CR protons r ga i(E) oc E -a5±ai 
fl6fl , and has a weaker energy dependence than r ICS does. At sufficiently low energy, 
then, r gal < T[ CS , and processes other than ICS cooling become relevant. The slope 
of the electron spectrum of Eq. (^) should change as the energy is lowered. The 
observed spectrum of Fig. 1 shows such a change, but it occurs below E e = 10 GeV, 
a range in which local modulations would mask the effect. 



Data from the SAS 2 satellite |L9| suggested the existence of an isotropic, 'dif- 
fuse' gamma background radiation (GBR). The EGRET instrument on the Compton 
Gamma Ray Observatory has confirmed its existence. By removal of the Galactic- 
disk diffuse emission a uniformly distributed GBR has been found, of alleged ex- 
tragalactic origin |J. The GBR flux in the observed energy range of 30 MeV to 
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120 GeV is well described by a single power law: 

-2.10±0.03 



dF 7 



f2.74 ± 0.11) x 10" 



MeV 



cm- 2 s" 1 sr" 1 MeV" 1 . (6) 



Many possible sources for the GBR have been discussed, and are reviewed in 0. 

EGRET data on 7-rays above 1 GeV show an excess over the expectation from 
cosmic-ray production of 7r°s [ 20| . Electron bremsstrahlung in gas is not the source of 
the 1-30 MeV inner-Galaxy 7-rays observed by COMPTEL [2lj|, since their latitude 
distribution is broader than that of the gas. These findings imply that a source such 



as ICS may be more important than previously believed [22]. We go even further 
and propose that the GBR is not mainly extragalactic, but dominantly due to the 
ICS of CBR and of starlight photons by CR electrons in the galactic halo. It is 
useful to discuss the CBR and starlight contributions in turn. 

Let us characterize the CR electron density in the Galaxy and its halo as a smooth 
spherical distribution. Assume the electron flux, Eq.(Q), observed at E e > 10 GeV, 
to be representative of the local interstellar value. Define a 'scale radius' for the CR 
electrons, R e = J dr[dn e (r)/dr]/n e (0). Let 9 be the angle away from the direction 
to the galactic center, distant d ~ 8.5 kpc. The CR-electron density distribution 
in energy and position may be convoluted with the CBR spectrum to derive the 
resulting flux of ICS up-scattered photons p3| . For R 2 3> dg the result is insensitive 
to the spatial shape of the electron distribution and very well approximated by: 



dF 7 n a T R e f(0) 



dE„ 



2 E. 
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E e = m e c \ I - — 
4e 



(7) 



where E e is obtained by inverting E 7 . The form-factor f(8) reflects our excentricity: 



R e f(0) = d cos 6 + \/Rl - d| sin 2 



(8) 



For the value we shall estimate, R e = 30 kpc, f(6>) = 1 ± 0.28, for 9 = 0, it. 



Substitute Eq. (g) into Eq. (J7p to obtain: 



dF 



CBR 



dE, 



1.43 ±0.30) xlO" 



30 kpc 



MeV 



-2.10±0.05 



cm-V^r^MeV -1 . (9) 



The index coincides with the measured one and, for R e = 30 kpc, the normaliza- 
tion is roughly half of what is observed, Eq.(|6|). The electron spectrum of Eq. (Q) 
describes the data in the range E e > 5 GeV, so that Eq. (||) should be valid above 
E 7 ~ 100 keV. We have predicted a steepening of the electron spectrum at E e ~ 1.6 
TeV. By the same token, the spectral index (3 <~ 2.1 of Eq. (||) should increase by a 
currently unobservable A/? ~ 1/8 at E 7 ~ 10 GeV. 



It is difficult to model in detail ICS on starlight [|4j] 0. We make a coarse 
estimate of the effect, at high latitudes, of an electron distribution with large R e . 
Approximate the Galaxy's starlight, of average energy ~ 1 eV, as that produced 
by a source at its center with the galactic luminosity L t = 2.4 x 10 10 L . Assume 
that R 2 >> d and replace in Eq. (||) the halo column density of CBR photons in 
the 6 direction, n o R e f(0), by the same quantity for starlight: 



47rce^ 



dx 



* JO u 



d% — 2xd cos# + x 2 47rce* 



vr- 



d sin 9 



(10) 



7 



to obtain 
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-i -2.10±0.05 



dE. 



(1.14 ±0.24) xlO 



vrsin^ 



MeV 



cnrW^MeV -1 . (11) 



The point source approximation is not accurate enough for directions close to the 
galactic disk and galactic bulge, but the GBR data [|| mask |b| < 10° as well as 
|b| < 30° for |1| < 40°. 

The angular average of the sum of CBR and starlight contributions to the GBR, 
Eqs.(f|) and (pTTj) , is shown in Fig. 2. It has the observed spectral shape and, for 
R e = 30 kpc, the observed magnitude. The fitted value R e ~ 30 kpc is imprecise: 
the starlight to CBR ratio is proportional to e*/eo raised to a poorly determined 
power, 0.10 ± 0.05. Our spectral index is independent of direction, as observed, see 
Figs. 5 and 6 of [[J. The spectral- index data are given in || for 36 (b, 1) domains, in 
Figs. 3-7 we have redrawn them as functions of 9: 9 values for each half hemisphere. 
The statistical test for a flat distribution is surprisingly good: x 2 /d-o.f. ~ 0.5. 

The intensity of the GBR diminishes as the angle away from the galactic center 
increases, see Figs. 5 and 6a of 0, which we have redrawn as functions of 9 in 
Figs. 8-12. Even after eliminating the four data points closest to the galactic center 
(at 9 = 49.2°), about which the observers feel unsure ||, xVd.o.f. = 1.8 for constant 
intensity. The sum of Eqs. (^) and ( |TTD fits better the observed angular trend of the 
GBR, as shown in Figs. 8-12 (x 2 /d.o.f. = 1.3). This agreement would presumably 
further improve with a more realistic treatment of starlight. 
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For R e = 30 kpc, the luminosity of our galaxy in 7-rays of energy above E is: 



L E ~ (1.5 ±0.5) x 10 



10 



g 1 -0.10±0.05 



MeV 



erg s" 1 . (12) 



This result is slowly converging and uncertain, but a future 7-ray telescope, such as 
GLAST, could possibly see the corresponding glow of Andromeda's halo. 

To estimate the contribution of external galaxies to the GBR, some concepts 
and numbers need be recalled. Hubble's 'constant' is H = 100 h km s -1 Mpc -1 , 
with h ~ 0.65; Q m and Q\ are matter and vacuum cosmic densities in critical units; 
Q = Q m + Q\; y = 1 + z is the redshift factor. In a Friedman model, the time 
to redshift relation is dy/dt = -H f(y) y, with f(y) = [(1 - Cl) y 2 + fi m y 3 + ^a] 1/2 - 
The luminosity density of the local universe is p L ~ 1.8 x 10 8 h L Q Mpc~ 3 . The 
combination p L /L^ provides an estimate of the average number density of 'Milky- 
Way-equivalent' galaxies. If the sources of CRs are young supernova remnants or 
gamma-ray bursts, the CR production rate ought to be proportional [^J to the star 
formation rate Rsfr, recently measured up to redshift z ~ 4.5 [f27 . 



The energy of CBR photons up-scattered by electrons at 'epoch y' is propor- 
tional to T(y) = y T and is subsequently redshifted by the same factor; hence the 
corresponding spectra from distant galaxies have the same energy dependence as 
from our galaxy up to E 7 ~ 10 GeV, where absorption on the infrared background 
becomes relevant f2"gfl. The cosmological starlight contribution is redshifted and 
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consequently becomes negligible. For the sum of all galaxies, we obtain: 

4ttR2 p L c f Rsfr(y) y dy" 



dF 7 ^ dFg R 
dE 7 ~ dE 7 



1 + 



(13) 



3 U H J x Rsfr(O) f(y) y 3 
where the CBR flux is that of Eq. (|9|). For Rsfr we interpolate the summary 



values of J27|. For Vt = fi m = 1 and R e = 30 kpc, the extragalactic contribution in 
Eq. ( |T3"D is at the 4% level (the h-dependences cancel). Other sensible values of the 
cosmological parameters give similar results. 

We have given a simple explanation for the observed spectral index of CR elec- 
trons. We have suggested that these electrons populate a "halo" region well beyond 
the disk of our Galaxy and that they produce the GBR by inverse Compton scatter- 
ing on CBR and starlight photons. Thus, the diffuse 7 radiation (in directions other 
than those of the Galaxy's disk and bulge) would not originate outside the Galaxy, 
but outside its starry domain. We have demonstrated that the predicted spectral 
index and angular distribution of the produced radiation agree with the GBR data, 
their normalization results in an estimate of the scale-height of the electron dis- 
tribution. Our results can be checked in various ways. The power index f3 of the 
electron spectrum should steepen above E c w 1.6 TeV by A/3 ~ 1/4. The halo of 
Andromeda should shine in gamma rays above a few MeV. The GBR should reflect 
the asymmetry of our off-center position in the galactic disk and its energy spectrum 
should not have the GZK-like cutoff expected [^] for cosmological sources. 
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Energy (MeV) 



Figure 1: The primary cosmic ray electron spectrum as measured || by Prince 
[crosses], Nishimura et al. [squares], Tang [circles], Golden et al. [triangles] and 
Barwick et al. [stars]. The theoretical curve is normalized to the data. 
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Figure 2: Comparison between the spectrum of the GBR, measured by EGRET H, 
and the prediction for ICS of starlight and the CMB by CR electrons. The slope is 
the central prediction, the normalization is the one obtained for R e = 30 kpc. 
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Figure 3: EGRET data for the angular dependence of the GBR spectral index as a 
function of direction, plotted as a function of 9: the angle away from the direction 
of the galactic center. The galactic disk and galactic bulge are masked: the data 
are for |b| > 10° for |1| > 40° and |b| > 30° for |1| < 40°. The four points at each 9 
correspond to the four half-hemispheres. The straight line is the predicted spectral 
index. 
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Figure 4: Same as Fig. 3 for the half hemisphere with b > 0, 1 > 0. 
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Figure 5: Same as Fig. 3 for the half hemisphere with b > 0, 1 < 0. 



18 




100 120 

6 (degrees) 



180 



Figure 6: Same as Fig. 3 for the half hemisphere with b < 0, 1 > 0. 
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Figure 7: Same as Fig. 3 for the half hemisphere with b < 0, 1 < 0. 
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Figure 8: EGRET data, masked and organized as in Fig. 3, for the dependence of 
the GBR intensity on 9. The line is the prediction for ICS of starlight and the CMB 
by CR electrons, for R e = 30 kpc. 
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Figure 9: Same as Fig. 8 for the half hemisphere with b > 0, 1 > 0. 
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Figure 10: Same as Fig. 3 for the half hemisphere with b > 0, 1 < 0. 



23 



-4 

xlO 




Figure 11: Same as Fig. 8 for the half hemisphere with b < 0, 1 > 0. 
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Figure 12: Same as Fig. 8 for the half hemisphere with b < 0, 1 < 0. 



25 



